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1.0 QHfiflDUGIIQN 

Three studies are presented in this document. The first is a followup of a pre- 
vious report (ref. 1), which presented a new technique for predicting monthly 
sunspot numbers over t period of several years (February 1979 to 1989). This 
followup report presents the result of a subsequent effort to further verify 
the prediction technique over the period for the maximum epoch of solar cycle 
19. The second and third studies presented in this report, respectively, are 
new techniques and results for predicting solar flux (F10.7 cm) based on two 
flux/sunspot number models, ascent and descent and geomagnetic activity indices 
A p and Kp as a function of sunspot number and solar cycle phase classes. 

Orbit lifetime models used within the MPAD generally utilize atmospheric density 
models that require data specifying a measure of solar flux and the geomagnetic 
activity index (fig. 1). 


2.0 SUNSPOT NUMBER PREDICTION TECHNIQUES 

The sunspot number prediction technique by the author and G. G. Johnson (ref. 
1), is applied here for additional verification. A different timeframe (years 
of the solar cycle) was used, which includes the maximum epoch sunspot numbers 
of solar cycle 19 (largest values observed to date). 


3.0 CYCLE. 19- SUNSPOT NUMBER PREDICTION 

The prediction analysis from (ref. 1) used the solar activity sunspot number his- 
torical period (monthly observed values) shown in figure 2 as input to the 
study. Also included in the document was a prediction verification analysis 
that used the 3ame sunspot number record but terminated the last input data 
point in December 1972, allowing actual observed sunspot numbers to be compared 
with the predicted results (January 1973 to January 1979). The timeframe 
included sunspot numbers observed just before and after the minimum epoch of 
cycle 21, June 1976. 

This :-'npariron showed highly favorable results for that time period, however, 
to fully test this prediction technique, it was deemed necessary to make this 
follow-on study using a different set of conditions; i.e., solar maximum epoch 
values. 

Presented in figure 3 are the observed sunspot numbers from approximately 1943 
to 1968, which includes the maximum epoch values this study will predict. The 
actual input data to the prediction model used are sunspot numbers observed from 
1830 to January 1956 (fig. 4). The resulting 10-year sunspot prediction starts 
in February 1955 and extends to January 1965. 


4.0 EIML-BESUBUAL CURVE 

By applying the same mathematical and computational technique of reference 1 to 
the sunspot number record in figure 4, a model composed of 19 spectural lines 
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was obtained by numerous iterations. Figure 5 shows the final residual ourve 
R19 as a function of month/years with the last value at January 1956. Using 
the same logic as before, the predicted residual was determined by symmetric ex- 
tension and is presented in figure 6. 

NOTE: The summation of the 19 harmonics of the model with the residual R1 r in 

figure 5 yields the original input sunspot record. Therefore, as pointed out in 
reference 1 , when the extrapolated portion of the residual R19 is summed in a 
like manner, the result is a prediction of the sunspot numbers for that 
timeframe . 


5.0 COMPARISON OF THE PREDICTED AND OBSERVED SUNSPOT NUMBERS 

The predicted sunspot numbers from January 1956 to January 1965 are shown in fig- 
ure 7 with the actual observed sunspot numbers for comparison purposes. The cus- 
tomary format of presenting these values (all positive) is presented in figure 8. 


6.0 SUMMARY OF CYCLE 19 SUNSPOT NUMBER PREDICTION 

To summarize the results of the R19 residual study, three plots are shown in 
figure 9 that include (a) the residual R 19 curve alone, (b) the R19 residual 
curve and extrapolated residual, and (c) sunspot number curve for both predicted 
and observed values. 

NOTE: These predicted sunspot numbers show good overall agreement with the ac- 
tual observed sunspot numbers. During the maximum epoch of this cycle, the maxi- 
mum delta sunspot number difference between the predicted and observed values 
was (±20) for a given month. However, when considering the entire curve for 
this timeframe, the predicted values average out to a respectable comparison. 
Therefore, the results from this R19 study indicate the predicted sunspot num- 
bers for the reference 1 study are acceptable. 


7.0 PREDICTED SOLAR FLUX (F10.7) 

The predicted solar flux data presented in reference 1 (1979-1989) were smooth 
F10.7-cm values based on an (sunspot number/flux) empirical relationship using 
smooth 13-month smooth sunspot numbers as input (fig. 10). The curve in figure 
10 is derived from a numerous amount of observed sunspot numbers and flux data 
without regard for the ascending or descending phases of the solar cycle; i.e., 
all of the data were included into one large sample. (See reference 1 of refer- 
ence 1 . ) 

This document presents the effects of the ascending and descending phases of the 
solar cycle and the two new (sunspot number/flux) models that were derived. 
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It was noted, when looking at the sunspot number/flux daily or monthly values 
over a period of time, that a distinction between the ascending and descending 
values existed and could be identified. The investigation of these data is 
discussed . 

a. Sunspot number and solar flux time history investigated 

Over 4000 observed daily sunspot number and flux data samples from 1955 
through 1965 (ref. 2) were used for the analysis. To understand the rela- 
tionship of sunspot numbers to flux values, several plots were generated and 
are discussed as follows: 

Figures 11 and 12 present, respectively, sunspot number and solar flux shows 
values observed over the 1 ,ear span. The same data (fig. 13) sunspot 
number as a function of flux values for this same timeframe. This plot 
shows all of the data samples together without regard to the ascending or 
descending effect of the phases of the solar cycles. The author 
investigated the difference between ascent and descent values. Presented in 
figure 14 are the ascent phase (sunspot number/flux) data extracted out of 
figure 13 and likewise the descent phase (sunspot number/flux) data shown in 
figure 15. The next step was to develop the ascent and descent models that 
were derived by determining an average flux value for a given sunspot number 
with its associated maximum and minimum values. The maximum and minimum 
values average out to be slightly less than ^2 sigma. This approach 
produced the models presented in figures 16 and 17, which are called ascent 
flux envelope and descent flux envelope, where the average value of each is 

included with its maximum and minimum expected values, thus determining the 

envelopes. 

b. Comparison of predicted and observed flux (1965-1979) 

A test case was generated (fig. 18) using the observed monthly mean sunspot 
numbers for (1965-1979) with the new models using predicted flux and 
observed flux for verification. The predicted flux agrees well with the 

observed flux and confirms the accuracy of the new models. 

c. Comparison of predicted flux data (1979-1989) 

This section presents various results of both the Johnson Space Center (JSC) 
and Marshall Space Flight Center (MFSC) predicted flux values for the years 
(1979-1989). The MSFC data used are from reference 3- Using the JSC 
predicted sunspot numbers from reference 1 and the new JSC au.ispot 
number/flux models (ascent and descent), a predirted flux envelope was 
developed and is presented in figure 17. Note that these flux values 
(maximum, average, and minimum) are predicted monthly mean values and not 
smooth as in the case for the MSFC prediction techniques. Figure 20 pres- 
ents the MSFC (5-month smooth) predicted flux envelope from reference 3 for 
cycle 21 where the maximum represents (97.7), minimum (2.3) and the average 
(50) percentile predictions, or approximately ±2 sigma values. 
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The JSC predicted flux me modified from the monthly values to a 5-month 
cuvve for comparison purposes with the MSFC 5-month prediotion. These data 
ave shown in figure 21 in conjunction with monthly observed flux. 

Finally, a summary plot for this flux study is presented in figure 22, whioh 
shews the JSC developed ascent and descent sunspot number/flux models and 
the predicted JSC monthly flux envelope for 1979-1989. 


8.0 GEOMAGNETIC INDICES Ap and K p 

Given a particular level of geomagnetic activity measured by the planetary geo- 
magnetic index Kp, the geomagnetic heating (increase in the exospheric tempera- 
ture above the quiet temperature corresponding to Kp * 0), can then be expressed 
as a function of the invarient magnetic latitude (ref. 4). 

a. Definition of A p and K p 

Excellent descriptions of these indices can be found in reference 5. Here 
it may be sufficient to say that the index Kp ranks all 3-hour intervals 
of the Greenwich day according to increasing disturbance (activity) by 
assigning one of the indices zero (quiet), 1, 2, 3, .... 9 (intense storm). 
The average of the eight ranges of (ap) per day is the daily amplitude Ap, 
which can range in value from zero to 400. The relationship between Kp 
and Ap . is presented in figure 23 for conversion purposes. 

b. Sunspot number/Ap prediction models 

Since the geomagnetic activity A p ex' K p indices predictions are used in 
orbital lifetime programs along with the predicted solar flux (F10.7 cm) 
values, the author performed a study to predict values of Ap corresponding 
tc> the already defined predicted sunspot numbers of cycle 21 (ref. 1). 

To derive suitable models for predicting values of A p , data were extracted 
from reference 6, which discussed and presented results of time variations 
of the geomagnetic activity indices Kp and Ap. Basically, this reference 
presented Kp and Ap data for the 12 calendar months using various 
classes in tne sunspot cycle for 30 years of data as follows: 

Class A - minimum epoch 
Class B - ascending 
Class C - maximum epoch 
Class D - descending 

For each class there exist*, a corresponding monthly Ap index that is the 
average for all days; the A p D index (average for 5 disturbed days per 
month and the ApQ index that is the average for 5 quiet days per month. To 
formulate a geomagnetic model for this study, the conoept of classes for the 
sunspot cycle was used similar to the above. However, the author redefined 
the classes for prediotion purposes as shown in figure 24. A description of 
these classes is presented in table I, and an example of one of the models 
(class B) is presented in figure 25. Each of the classes Las been modeled 
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as ApD, Ap and ApQ as a function of sunspot number as shown in this 
plot. Therefore, using the predicted sunspot numbers (monthly) from refer- 
ence 1, which includes solar cycle 21 and part of solar cycle 22, and these 
geomagnetic models (classes of the sunspot cycle), then predicted values of 
ApD, Ap, and ApQ were obtained for this same timeframe. 

c. Comparison of predicted geomagnetic indices 

The predicted JSC geomagnetic envelope values of ApD, A p and ApQ from 
1979 to 1989 are presented in figure 26. 

The predicted MSFC geomagnetic envelope values of A p from 1979 to 1987 from 
reference 2 are pi esented in figure 27. They are 5-month smooth Ap index 
values where the maximum represents (97.7), minimum (2.3) and average 50- 
percentile predictions, approximately ±2 sigma. 

Because the MSFC A- data are based on 5-month smooth values, an additional 
step was preformed to alter the JSC monthly Ap values to predicted 5-month 
smooth values. Comparison of the JSC 5-month smooth Ap values with the 
MSFC 5-month smooth envelope data are presented in figure 28. The JSC data 
are not as smooth as the MSFC but overall agree favorably. Also included is 
the JSC 5-month smooth Ap curve imposed on the plot in figure 26. 

In summary of these comparisons figure 30 is presented, which shows three sepa- 
rate plot cases. Included are the predicted MSFC 5-month smooth Ap envelope, 
the JSC monthly Ap envelope, and the same JSC envelope with the imposed MSFC 
5-month smooth Ap values plotted as connecting circles. 

Because some organizations use K p instead of A p as input to their density 
models, then converting to predicted Kp values is accomplished with little dif- 
ficulty using the relationship of A p to Kp already presented in figure 23- 
Figure 31 presents a plot of the predicted Kp values that were converted from 
Ap. The data shown are for the MSFC 5-month Kp envelope in conjunction with 
the JSC 5-month Kp values. 


9.0 CONCLUSIONS 

This paper deals with three primary subjects: (1) predicted sunspot numbers, 

(2) predicted solar flux (F10.7 cm), and (3) predicted geomagnetic indices Ap 
and Kp. The techniques and prediction models for each of these are presented 
and discussed with numerous plot results. 

To best summarize this work, included is a final summary plot (fig. 32), which 
presents the final predicted results of each of these studies. 
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TABLE I.- SUBDIVISION OF CLASSES FOR A TYPICAL SOLAR CYCLE 

Claas Description 

A Immediately before minimum epoch (aunspot number range ia 

approximately: 5-12) 

A' Immediately after minimum epoch (sunspot number range is ap- 
proximately: 2 -12) 

A+ After class A* (to sunspot number approximately: 35) 

B Ascending (sunspot number range to approximately: 93) 

C Ascending maximum epoch (sunspot number range to approxi- 
mately: 100) 

P Descending maximum epoch (sunspot number range to approxi- 
mately: 100) 

E Descending (sunspot number range to approximately: 26) 

E' Descending before class A (to sunspot number approximately: 

25.9) 
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Figure 1.- Flow chart for predicted solar activity Input values. 















SUNSPOT NUMBER CYCLES (18 RND 19) 



Figure 3.- Sunspot number time history of values to predict (1956-1965). 







Figure 4.- Sunspot number record used as Input to study (1830-January 1956). 
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Figure 5.- Final residual curve R19 (1830-January 1956). 
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Figure 6.- Final residual and predicted residual curve R19. 
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Figure 7.- Sunspot numbers (observed and predicted) (1 943-January 1965). 








RESIDUAL 




0009- 

oooe- 

0 

oooe 

0009 

tn 

o **« 

O 

8 

I 

9000 

Q O O O 
Q Q O O 
O O o ° ° 
» n fn «) 

t i 

cn 

O w 
O 
O 
Ok 

t 

oo o io 

iM ru (\i h 

f\J o o ° 

H 00 t 




ckrujco^Qzxr^-J 


irorincLo*- 

ZDCMUtt 


NTH'YEARCS) 

OBSERVED AMD PREDICTED SUNSPOT NUMBERS 













SUNSPOT NUMBER/SOLflR 


} 



wo z'ou xnu yyios nuow hioows 


17 


SMOOTH MERN ZURICH SUNSPOT NUMBER 






<400 



18 


1010 2020 3030 4040 

DAILY OBSERVATIONS 11 YEARS 











Figure 13.- Solar flux as a function of sunspot numbers for all samples. 
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function of sunspot numbers. 






22 



PLUX/SSN RSC . ( ENV . ! - JSC MODEL 





I 

a 

i 

o 

i 

o 

i 

a 

i 

o 

i 

o 

1 

o 

o 

• 

s 

d 

in 

• 

o 

o 

d 

in 

d 

o 

• 

o 

LO 

d 

o 

s 


to 

ro 

CM 

ro — 





samuA'z 

•ou'xmj 




Figure 16.- ascent flux/sunspot number envelope model 
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Figure 17.- Jsc descent flux/sunspot number envelope model. 
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Figure 19.- JSC predicted monthly solar flux envelope (1979-1989) 
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Figure 20.- MSFC predicted 5-month smooth solar flux envelope (1979-1987). 
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Figure 22.- Predicted solar flux envelope with the 
JSC sunspot number/flux models. 
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Mgure 24.- Classification of the solar cycle for the geomagnetic index A D models. 
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Figure 26.- JSC predicted envelope of the mean geomagnetic 
indices (A d D, A d , A d Q) for (1979-1989). 
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Figure 27.- MSFC preoicted envelope of the 5-month smooth 
geomagnetic index A_ (ly79-1987). 
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Figure 29.- Comparison of the MSFC 5-month smooth Ap predict 
with the JSC predicted monthly A D envelope. 
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Figure 30.- Summary plots of the MSFC and JSC Ap envelopes. 
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Figure 31.- Predicted MSFC 5-month K p envelope and JSC 5-month 
average values (1979-1939). 
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Figure 32.- Summary plots of predicted sunspot numbers, 
predicted solar flux (F10.7 cm), and predicted A n , );<r 
values (1979-1989). ” f ' l/ C& V 
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